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Abstract 

We consider the leading non-analytic temperature dependence of the specific heat and temper- 
ature and momentum dependence of the spin susceptibility for two dimensional fermionic systems 



'^ . with non-circular Fermi surfaces. We demonstrate the crucial role played by Fermi surface curva- 

O I ture. For a Fermi surface with inflection points, we demonstrate that thermal corrections to the 



uniform susceptibility in D = 2 change from Xs oc T to Xs oc T"^'^ for generic inflection points, and 
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^ ■ to Xs oc T^' ^ for special inflection points along symmetry directions. Errors in previous work are 

\o _ 

0> , corrected. Application of the results to Sr2Ru04 is given. 
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I. INTRODUCTION 

L. D. Landau's "Fermi liquid theory" provides a robust and accurate description of the 
leading low temperature, long wavelength behavior of a wide range of systems of interacting 
fermions in two and three spatial dimensions. In Landau's original worki it was assumed 
that the temperature (T) and momentum (g) corrections to the leading Fermi liquid behavior 
were analytic functions of (T/Tp)^ and [q/kpY, with the Fermi momentum kp set by the 
inter-particle spacing and the Fermi temperature Tp ~ vpkp with vp a typical measured 
electron velocity. However, subsequent work revealed that in dimensions d = 2 and d = 3 
the leading temperature and momentum dependences of measurable quantities including the 
specific heat and spin susceptibility are in fact non- analytic functions of T^ and g^. The 
nonanalyticities have been studied in detail for Galilean-invariant systems (spherical Fermi 
surface or circular Fermi line) (see Ref y for a list of references). In this paper we extend 
the analysis to a more general class of systems, still described by Fermi liquid theory but 
with a Fermi surface of arbitrary shape. 

The extension is of interest in order to allow comparisons to systems such as Sr2RuO^ 
and quasi-one-dimensional organic conductors^, in which lattice effects are important. The 
extension also provides further insight into a fundamental theoretical issue: a crucial finding 
of the previous analysis^'^ was that for two dimensional systems the non-analytical term 
in the specific heat coefficient 6C{T)/T oc T arise solely from " backscattering" processes 
at any strength of fermion-fermion interaction. For the spin susceptibility, the situation is 
more complex: in addition to the backscattering contributions to 6xs(T) oc T and Sxsil) oc 
q there are contributions of third and higher order in the interaction which involve non- 
backscattering processes^^. If the interaction is not too strong the backscattering terms 
dominate. 

Unlike scattering at a general angle, the kinematics of backscattering is effectively one- 
dimensional, and depends sensitively on the shape of the Fermi surface. We shall show that 
the crucial parameter is the Fermi surface curvature, and that for two dimensional systems 
in which the Fermi surface possesses an inflection point, the power laws change from |(T, q)\ 
to |(T, g)^/^| or |(T, g)^/'^| according to whether the inflection point is or is not along a 
reflection symmetry axis of the material. Similar effects occur in three dimensional systems, 
but the effects will be weaker as the singularities there are only logarithmic. 




FIG. 1: Fermi surface of Sr2Ru04 computed from tight binding parameters deduced from quantum 
oscillation measurements^ , with a few "parallel tangents" points indicated by letters A,B,C,D and 
the Brillouin zone points M = (vr, 0) and X = (vr, vr) also noted. 

The importance of the Fermi surface geometry was previously noted by Fratini and 
Guinea, who showed that the presence of inflection points changes the power law for the 
spin susceptibility Xs{T)^- However, we believe that their calculation treated the kinematics 
of the backscattering incorrectly; as a result, they found that, in 2D, the presence of the 
inflection point changes the temperature dependence only by a logarithm, from Xs{T) oc T 
to XsiT) oc TlogT, instead of the Tv^'sj found here. 

The paper is organized as follows. Section II introduces the model and defines notation. 
In Section III, we demonstrate the physical origin of the results, via a calculation of the long 
range dynamical correlations of a Fermi liquid. Section IV presents results for the specific 
heat of a multiband system. In Section V we calculate the nonanalyticities in the momentum 
and temperature dependence of x- Section VI shows how new power laws emerge for Fermi 
surfaces with inflection points, and why one-dimensionality affects the powers. Section 
VII presents estimates of the size of the nonanalytic terms in Sr2RuOi. Section VIII is a 
conclusion. 



II. MODEL 

We study fermions moving in two dimensions in a periodic potential. As discussed at 
length in Ref. [2], because we are concerned with the low T properties of a Fermi liquid we 
may adopt a quasiparticle picture. Lifetime effects are not important, and the quasiparticle 
weight (z) factors may be absorbed in interaction constants. We may therefore consider 
several bands, labeled by band index a, of quasiparticles moving with (renormalized) dis- 
persion Ep. The Fermi surfaces are defined by the condition e^ = /i; an example is shown in 
Figd We parameterize the position at the Fermi surface by a coordinate s". For vectors k 
near a particular Fermi surface point k^ we will write 

eJ-, = 4M(.„ + ^) (1) 

with vf{s"') = \del_/d]<.\ the Fermi velocity at the point (s°) and the components of k parallel 
and perpendicular to the Fermi velocity (Fermi surface normal) given by 



A;,| = (k-k^)-9^(s'^) (2) 

k^ = {k- kp) X v^(s") (3) 

kQ^^s"") is the curvature of the Fermi surface at the point s°. For a circular Fermi surface, 
ko = kp independent of s; but in general ko ^ kp and both depend on s. 

It is sometimes convenient to use the variables e^ and ^fc, where Q\^ is the angle determining 
the direction of the Fermi velocity vrik) = dsk/dk relative to some fixed axis. Eq^ shows 
that the Jacobean of the transformation is 

d'k = ^^deudOu. (4) 

In a non-Galilean-invariant system the Fermi surface may contain inflection points at 
which the curvature vanishes, i.e. k^ -^ oo. If the inflection point does not lie on an axis of 
reflection symmetry of the Brillouin zone, the dispersion (measured in terms of the difference 
of the momentum from an inflection point) is 

Sk - Ai = ^^F ( /^ll + Tj ) (5) 



However, if the inflection point lies on a symmetry axis, then only even powers in k± may 
occur and 

Ek- fi = VF ik\\ + -^j (6) 

Here ki and k2 are coefficients expected in general to be ~ kp. 

The fermions interact. We assume that the T — i> 0, long wavelength properties are 
described by the Fermi liquid theory, so that at low energies the interactions may be 
parameterized by the fully reducible Fermi surface to Fermi surface scattering amplitude 
^a,f3,'ys{k,P', k,p). For particles near the Fermi surface, |k|, |p| ^ kp, and F depends on the 
angle 6 between k and p and on the band indices a (for k) and b (for p). Backscatter- 
ing corresponds to 6* = tt. It is often useful to decompose F^^^.^ 5(7?) into charge and spin 
components 

ra,/3;7,5(7r) = F" '""SajS/JS + F" 'Va^CT/j^ (7) 

It is also instructive to make contact with second order perturbation theory for a model 
in which the particles are subject to a spin- independent interaction 

i:^^n< = $^t/(q)PqP-q (8) 

q 

with charge density operator 



pa 



The leading perturbative result is then 



r = c/(o) - ^^PM. r- = -^^(HM (10) 

If the interaction is local and only one band is relevant, i.e., U{q) = U, V^ = — F* = U/2, 
and 

ra,/3;7,5(71") = U {6a-ySf3S - SasSjS-y) (ll) 

III. PHYSICAL ORIGIN OF NONANALYTICITIES; ROLE OF CURVATURE 

Previous studies of the isotropic case demonstrated^*^ that the nonanalyticities in the 
specific heat and spin susceptibility arise from the long spatial range dynamical correlations 
characteristic of Fermi liquids. These are of two types. One involves slow (|fi| < vpq) long 
wavelength fluctuations and is expressed mathematically in terms of the long wavelength 



limit of the polarizibility ^Ilivy = limg^o^iQ,^) — n(g, 0) ~ |^| /o'- The 1/q behavior 
of polarizibihty gives rise to a long-range correlation between fermions which decays as 
|fi|/r at distances 1 <^ rkp < Ef/\uj\. The other involves processes with momentum 
transfer q ^ 2pp, and is expressed mathematically in terms of the polarizibility 611 = 



n(g,fi) - U{2kF,n) ~ n/^2kF - q (for \Vt\ < 2kF - q. The l/V^kp - q behavior of 
611 gives rise to an oscillation with a slowly decaying envelope, cos{2kFr — ^T/4:)\Q\/^/r , 
again at distances 1 -C rpF < Ef/\uj\ and again leading to singular behavior. We now 
compute these processes in the multiband model defined above, and then show how they 
affect thermodynamic variables. 

Consider the long wavelength process first. The non-analyticity comes from a particle- 
hole pair excitation in which both particle and hole are in the same band, and have momenta 
in the vicinity of Fermi surface points s* satisfying vf{sD ■ g = 0. Choosing one of these 
points as the origin of coordinates we have 

dk\\dkj_ 1 



611' 



H ^U,,{q,n)-U,,{q,0) = Tj2f 



(2-)^ ^cOn-VF(s:)(h + ^^%;^ 



^ \ aJ \ II 2fco(s*) 



"" ' .-^ + ^^n - VFisl) {k, + i^0f) ..„ - VFis:) [k, + ^^0f) ' ^''^ 

Performing the integral over k\\ and the Matsubara sum as usual yields 

^n:? = / ^^- %-- (13) 

/ 27r 2nvFis*) iO _ "fJO'^^'J ^ ' 

We see that for the part of 611 which is even in Vt the integral is indeed dominated by 
fcj_ ~ r2fco(s*)/fF('5*), so the approximation of expanding near this point is justified, and 
we obtain the nonanalytic long wavelength contribution as a sum over all Fermi points s* 
satisfying vf{s*^ ■ g = with coefficients determined by the local Fermi velocity and local 
curvature: 

For a circular Fermi surface, two points satisfy vf{s*^) ■(] = 0, k^ = kp and Eq. (fT^ reduces 
to the familiar result kF\^\/{27TVpq). If the curvature vanishes, then use of Eqs (j3j) or (jU)) in 
Eq (fT^ yields a 6I[lw ~ i^/l) ' respectively. 

We next consider the "2fcir" process. Here the situation is a little different. The singu- 
larities in general come from processes connecting two Fermi points with parallel tangents 



(the importance of parallel tangents points has been noted in other contexts^)-for example 
the points A, B or A,C shown in Fig ^ For a given vector q we denote as Q the closest 
vector which is parallel to q and connects two "parallel tangents" points. Symmetry ensures 
that the leading dependence of ^IIq = n(q) — n(Q) involves only the gy = (q — Q) ■ Q/|Q|. 
Labeling the initial and final of the two Fermi points connected by Q as Si 2 and noting that 
for systems with inversion symmetry the points come in pairs symmetric under interchange 
to the band indices we have 



f, . X /„ X v^ f dk\\dk\ 1 

OT J ^ n(q, n) - n(Q, 0) = T 5^ / ^iL^- --^ 



iuJn -Vf{Si) (/C|| + 
1 1 



iUn + i^n - Vf{s2) U|| + q\\ + 2ka(s2) ) ^^" ~ Vf{s2) [k\\ + 




(1^2115) 



Here the sign of vp and k^ becomes important. For two points "on the same side" of the 
Fermi surface (e.g. points A and B in FigQ) the two velocities and the two curvatures have 
the same sign (a change of k\\ either increases or decreases both energies), the integrations 
proceed as in the analysis of Eq (fT^ . and the different position of the q means that to the 
order of interest there is no singular non-analytical term. On the other hand, if the two 
velocities have opposite sign (e.g. points A, C in Fig^ then after integrating over fcy and 
A;_L we obtain (for fi > 0) 

5Ii^= ,f^\ T y _J£±^ + (1^2) (16) 

Here vf\ = vf{si), vf2 = vf{s2), and 

1 1 



kavg 2 \A;o(si) kQ{s2 

1 1/1 1 



Vavg 2 VI^Fll \vF2 

Completing the evaluation yields 



+ 77^ ) (17) 

(18) 




The singular |r2|/y/g behavior of the dynamic 6IIq only holds when gy < and is obtained 
by expanding (fT^ in ^/q\\ ■ On the contrary, the singular behavior of the static Hg oc ^/gj[ 
holds at gy > 0. The dynamic part of AHq behaves as f2^/(g||)'^/^ at gy > 0. 




FIG. 2: One of the two second-order diagrams for the free energy, which give the non- analytic 
contribution to the specific heat (from Ref:2j). Here T is the (unsymmetrized) fuhy renormaUzed 
Fermi surface to Fermi surface scattering amplitude. 



IV. SPECIFIC HEAT 



This section treats the non-analyticity in the specific heat. The GaUlean- invariant case 
studied previously is simple enough that the corrections can be evaluated, with no approxi- 
mations beyond the usual ones of Fermi Liquid theory^. The evaluation confirms that in two 
dimensions the non-analytical contributions to the specific heat involve only the backscat- 
tering amplitude r(7r). In work prior to that reported in Ref |2] this conclusion was reached 
by approximate calculations in which it was assumed that the non-analytical contributions 
were governed by the backscattering amplitude only, and then the assumption was shown 
a fortiori to be consistent. In the non- Galilean-invariant case of interest here, a complete 
analysis along the lines of Ref [2] is not possible. We will follow earlier work and assume 
that the effects arise only from backscattering processes, and then show that the assumption 
is self consistent. We specialize for ease of writing to a momentum-independent vertex F 
(see (jllj) ). but keep band indices. For a two-dimensional system the diagram which gives 
the nonanalytic term in the thermodynamic potential S is^ shown in Fig |21 We may write 
the resulting diagram schematically as 



IJ2 fidqdn)Tl^Ua,iq,n)I[Uq,^) 



abed 



(20) 



Now, singularities leading to non- analytical terms may arise for g ^ 0, in which case we 
must consider only the intraband contribution to 11, and q —>■ Q, where Q is one of the 
"parallel tangents" vectors mentioned in the previous section, in which case the band indices 
of polarizabilities may be different. Let us consider first the small q singularities. We have 



-^ E J {dqdn)Tl,,Ua,{q, n)Il,,{q, fi) (21) 

Substituting from Eq (fT^ gives 



■^LW 

ab 



^LW 



2 ^J (2vr)2Z^ ^"'^ q' {^^?viy,, ^ ^ 



The integral over q is logarithmic and is cut by Q; the analytical continuation and integral 
of frequencies may then be performed and we obtain 

3C(3)^ fderl,Me)k',{e) 



6C_ 



LW ""^ ^-^ 27r v'p,^Vp^^ 






We now turn to the parallel tangents part of the calculation, finding 

-Q = -IY1 j {dqdn)Vl,,JiM + q. ^)^UQ + g, fi) (24) 

ab 

Note that the symmetry of the vertex means that F^^^^ = Vj^^bb- Again substituting AIIq 
instead of Ilab{Q + q-, ^) and evaluating the integrals explicitly we find ® (note that to obtain 
the non- analytical behavior it is sufficient to expand 5xq for \VL\ « q) 

^ __}_\^r^ rr\^ f dq\\dqL ^^ kg^g , . 

-Q- 2 2_. ^ ^^^a^ Z. J (2^)2 (Anyv^^y,, \q, | ^''^ 

Now noting that kavg = ^0^0/(^0 + ^0) ^^^ ^^^^ dq± = (/cg + kQ)d9 we see that Eq and E^w 
give identical contributions despite the apparently different kinematics. 
Adding the two contributions, we obtain 

6C _ 3C(3) ^ f de Tl,Me)k',{e)_ ^^e) 



T 27r3 ^y 2n v\,v\. 

Note that the integrals over the Fermi surface contain k\ = kQ{6k) rather than the product 
of two fco factors at different points along the Fermi surface. This is a direct consequence 
of the fact that only backscattering contributes to (jinjl and (jHj). For a one band model 
with an isotropic Fermi surface, Eq. ()26|) reduces to the result in^. For a generic interac- 
tion, the calculation goes through as before with r^^j,^ replaced the components of the fully 



renorinalized, symmetrized Fermi surface to Fermi surface backscattering amplitude so that 
5C _ 3C(3) ^ f de (r«''-'^(7r)2 + 3r"''--(7r)) ^fcg(g)fcg(g) ^^^^ 



ab 



T 27r3 ^j 271 vlv% 



F,a^Ffi 



V. SUSCEPTIBILITY 
A. Overview 

This section presents calculations of the non-analytical momentum and temperature de- 
pendence of the spin susceptibility of a two dimensional non-Galilean-invariant Fermi liquid 
system, with a Fermi surface without inflection points. In contrast to the specific heat, there 
are two classes of contributions to the nonalnayticities in the susceptibility^^. One is of sec- 
ond order in the fulUy renormalized interaction amplitudes, involves only the backscattering, 
and is treated here. The other, which we do not study here, is of third and higher orders, 
and involve averages of the interaction function over a wide range of angles (analogously to 
similar contributions to the specific heat of a three dimensional Fermi liquid^). The former 
process is dominant at weak coupling, and involves an integral of the square of the curvature 
over the Fermi arc. The latter process has a less singular dependence on the curvature. 

Even to the order at which we work, many diagrams contribute (see FigE)); we evaluate 
one in detail to illustrate the basic ideas behind the calculation and then simply present the 
result for the sum of all diagrams. Consider for definiteness the "vertex correction" diagram 
- diagram 3 in Figure El The analytical expression corresponding to this diagram is (we use 
a condensed notation in which (dk) stands for an integral over momentum, normalized by 
(27r)^, and a sum over the corresponding Matsubara frequency) 

5x{q, 0) = -4 ^ [{dk){dl)G\k + q)G\k)Kf{l)G\k + / + q)G\k + /) (28) 

ab -^ 

and A is the product of the interaction vertices and internal polarization bubble: 

Af (g) = j{dp) {T'^'''\e)fG\p + q/2)G\p-q/2)] (29) 

where 9 is the angle between k and p, and we have used the fact that k and p are near the 
Fermi surface. A complete calculation in the Galilean-invariant case shows that, just as for 
the specific heat, the non-analytical momentum dependence of the susceptibility arises from 
the regions of small / and / ~ 2kp. We consider these in turn. 

10 








5. 





FIG. 3: Relevant second-order diagrams for the spin and charge susceptibihties (ferm Refj^). The 
last two diagrams are non-zero only for the charge susceptibility. 
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B. Small 1 

Here we choose a particular point Sk on the Fermi surface and integrate over Sk and the 
corresponding frequency. We parametrize the position on the Fermi surface by the angle </> 
between v^ and q and use Eq|l| We adopt coordinates /y and l± denoting directions parallel 
and perpendicular to v(sfe) and obtain 

1 1 '^ \ Ci \ 

7r-VA;-q-/|| -s- + Vfc-q-i|| -^ - Yk ■ q- l\\ I 

'k " "k " "k " / 

We may similarly evaluate A, proceeding from Eq (j29|) . Choosing as origin the point p = — k, 
defining coordinates p\\ and p_L antiparallel and perpendicular to v^, integrating over p\\ and 
the corresponding loop frequency gives 

Af(^ih^x,fi) = hr:r^ I ^ .. , r .,. (31) 



Yaabb^^Y 


f dp± iQ 


MvlY J 


1 O-rr i^ 1 7.. P±l± 
^k '^0 



Viewed as a function of l\\ the second hne in Eq flHUjl decays rapidly (~ IT) at large l\\ and 
has poles only in the half plane sgnlml\\ = sgnQ. Evaluation of the l\\ integral by contour 
methods, closing the contour in the half plane sgnl\\ = —sgnQ shows that nonanalyticities 
can only arise from singularities of A. Reference to Eq (J3T| shows that these can only arise 
from momenta p satisfying Vp ■ v/j < 0. In the Galilean-invariant case the integral could be 
evaluated exactly; the resulting non-analytical terms were found to be determined by a very 
small region around p = — k, i.e., around ^ = vr in ()29p. Here we assume that this is the 
case, and show a fortiori that the assumption is consistent. 

Performing the integral over Zy yields 



(32) 



2^ {2.f J 2n{vpir ^, 

\ ^avg Kq Vavg 


r dl^dpx_ ( ^ \ 
1 (2vr)2 Wk-q) 

2 


irp± 2iQ i±pj_ 



with Vavg defined in Eq^J Performing the sum over frequency and rescaling each of l±,p± 
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by y /cq yields 



^3 
^avg 

\2 (ya\2 



^^^^^^ + ... (33) 



ab 

'J{l±,P±; V ■ q) + r{h,pr, -V ■ q) - 2T(/^,p^; 0)) 



IVfc ■ qj 

with 

T(x, y; 2;) = (xy — z)'^ log |xy — z\ (34) 

Eq (|33|) is based on an expansion for small 1_l,p_l. The integrals over these quantities are 
cut off by other physics above a cutoff scale A which we have written as a hard cutoff. The 
ellipsis denotes other terms arising from physics at and beyond the cutoff scale, which lead 
to additional, regular contributions to 6x involving positive, even powers of q. Evaluation 
of Eq ESI yields (details are given in Appendix \K\i 

^ (r-^^(7r))^|q| f dd.k^X ^la 1^ ^^ r-,.^ 

SXLw{q) = -y- T^ / o / a\2f b\2f a\2 Vfc-? +-• (35) 

where Vk = vpiOk), k^ = ko{Ok), v^-g = cos(^fe — 6'^), and 9q is the angle between the direction 
of q and the direction of ^ = 0; the ellipsis again denotes analytical terms. We see that the 
non-analytical term is explicitly independent of the cutoff, confirming the consistency of our 
analysis. An alternative evaluation of Sxiwil) is presented in the Appendix |Bl 

For a circular Fermi surface, ko = kp, vp is a constant, and Eq. (jHSI) reduces to the 
previously known resullj^. However, the previously published computations are arranged in 
a way which apparently does not invoke the curvature at all. In Appendix C we show that 
the previous method does in fact involve the curvature, and leads to results equivalent to 
those presented here. 

C. 2pp processes 

To evaluate the contribution of 2pF processes we could proceed from Eq (j^ but ex- 
panding A in q = q — Q, where, we remind, Q = 2kFq- As before the products of G 
produce an expression with all poles in the same half plane. Exploiting the non-analyticity 
of n(q + q, fi) we obtain an expression which has a nonanalytic part which evaluates to the 
same expression as Eq (jSHl)- Instead of presenting the details of this calculation, we present 
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an alternative approach due to Belitz, Kirkpatrick, and Vojta®, in which one partitions the 
diagram into two " triads" , using the exphcit form of A^ 

Sx{q) = -4J2 f{dl){dki){dk2)T'"'''\ef[G\ki + q)G\k{)G\ki + l)] [G\k2 + q)G\k2)G\k2 + I)] 

ab -^ 

(36) 
where 9 is the angle between ki and k2. Choosing a particular point on the Fermi surface 
and evaluating the integral over e^i and the associated frequency yields {q\\ is the component 
of q parallel to the direction chosen for ki) 

ko{ei)dei Q 



Sx 



2Tr{v-y (v-.q) \^(f -/||) f -v-q-/,, 

^aabb^0^^2 ^Qb^j^^ ^ q)G\k2)G\k2 + /)] (37) 



As in the previous calculation, singular contributions can only come from regions where k2 
is directed oppositely to ki. Choosing as origin of /c2 the point diametrically opposite /ci, 
introducing parallel and perpendicular components as before and integrating over k2\\-, the 
associated frequency, and l\\ we get 

k'^{ei)ddi T\7t) ^^ f dk^dl^ f fi2ggj^^ 



^ 2^ J 27r(t;«)2(t;^)2 (27r)2 2^ J 



1 2 



«f^ — -Cr <-. _ ^xfc± 2in I - _, _ l±k± 2iO. Ilk, I \^°) 

„ V ■ q i, „ "T ^ ■ q i.b ., 

Eq ()38p is seen to be of precisely the same form as Eq ()32|) and gives the same result; the only 
difference is the dependence on orbital index. Integrating over frequency and combining the 
results from small q and 2A;^ contributions gives an answer whose orbital dependent part 
depends on the velocities via the combination 

v^iv^f + {v'^fv'' + 2(t;")2(t;^)2 



^ya _^ yb^ 



(v ■ q) = Vavg^ ■ q (39) 



D. Final Result 



Collecting the small q and 2kF contributions from all diagrams in Fig. El we find for the 
spin susceptibility 



ab 
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At g = and T > 0, we have 

^^^(^^ = 1^-3 I 2n{v^nvr ^^ (")))• (^^) 

For an isotropic Fermi surface and one band, this again reduces to the result in^. To make 
contact with previous work, which considered a simphfied interaction with indentical spin 
and charge components, r"'''*(7r) has to be replaced by r"''(7r)/2. Higher order powers of F 
do contribute to |g| and T terms in Xs, in distinction to C(T)/T, and at strong coupling, 
the non-analytic terms in the spin susceptibility are not expressed entirely via F^(7r)^. 
the terms of order F'^ have a less singular dependence on the curvature. At weak and 
moderate coupling, though, Eqs. ()40|) and ()41|) should be sufficient. Finally, for the charge 
susceptibility Xc{(l,T) non-analytic contributions from individual diagrams are all cancelled 
out: the full Xc{(l,T) is an analytic function of both arguments. 

VI. FERMI SURFACES WITH INFLECTION POINTS 

We see from Eq. fl40|) and (^1]) that as long as ko is finite all along the Fermi surface, the 
anisotropy of the Fermi surface affects the prefactors for |g| and T terms, but do not change 
the functional forms of the non-analytic terms in the specific heat and spin susceptibility. 
New physics, however, emerges when the Fermi surface develops infiection points at which 
ko{dk) diverges. Infiection points are a generic feature of realistic Fermi surfaces of two 
dimensional materials. In this section we show how inflection points emerge and then indicate 
the modifications they make to the results presented above. 

A. Inflection points in commonly occurring models 

We first note that many quasi-one dimensional organic materials have a band dispersion 
described by 

Horganic = -2ta {\kx\ - kp) - 2t'cos{kyb) (42) 

with |t'| << t, lattice constants a,b not too different, and a third dimensional coupling 
weaker than t' by an order of magnitude. In this case ko = 2t'cos{kyb) + 0{t''^ /t) obvi- 
ously vanishes at kyh k, ±7r/2. Thus inflection points are generic to quasi one dimensional 
materials. 

15 



We now consider the fully two dimensional t — t' model, with quasiparticle dispersion 

Ek = — 2t(cos kx + cos ky) + 4t' cos k^ cos ky — fi (43) 

We assume that t and t' are positive and /i is negative; t should be larger than 2t' for 
stability. 

Consider now the Fermi surface crossing along the (0, 0) -^ (vr, 0) direction (if it exists). 
The crossing occurs at 

-(2t - At')cos{kx) = fi + 2t (44) 

the velocity is along x and the curvature may be read off by expanding to second order in 
ky, giving 

^°(^) = {2t - Jcosik^)) ('') 

which is manifestly positive. 

On the other hand, at the Fermi surface crossing along the diagonal kx = ky we find 

which is positive at small k^. but changes sign as the Fermi surface approaches the point 
(7r/2,7r/2). We therefore conclude that for chemical potentials in the appropriate range 
inflection points must exist because the curvature has opposite sign at two points on the 
Fermi surface. 

B. spin susceptibility and specific heat 

We now analyze how the inflection points affect the nonanalytic terms in the spin suscep- 
tibility and specific heat. For simplicity, we restrict to one-band syatems. Quite generally, 
near each of the infiection points ko{6) behaves as 

ko{9) oc (^ - ^o)"' (47) 

At ^ = ^0) the curvature diverges, i.e., there is no quadratic term in the expansion of the 
quasiparticle energy in deviations from the Fermi surface. In the generic case, the dispersion 
is then 

£k = VFiOo)kii + Akl (48) 
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where, as before, the directions k\\ and k± are along and transverse to the direction of the 
Fermi velocity at ^ = ^o- For a special situation when 6*0 coincides with a reflection symmetry 
axis for e^, i.e., when vpi^^o) is directed along the Brillouin zone diagonal in t — t' dispersion, 
the expansion of Sk in the direction transverse to the zone diagonal holds in even powers of 
k±, i.e., at ^ = ^o; 

Sk = VF{eo)k\\ + Bk\ (49) 

In both cases, a formal integration over 9 in Eqs. (pn|) and ()41|) yields divergences. The 
divergences are indeed artificial and are cut by either A oi B terms in the dispersion. The 
effect on 6x{(l,T) and 6C{T)/T can be easily estimated if we note that the angle integrals 
diverge as j dOk^iO) ~ f d9/{9 — 6qY. In a generic case, described by (jlSj) . l/\6 — 6q\ has 
to be replaced by l/[\6 — 6q\ + Ak^^ty^. The angle integral then yields l/k±^typ. It follows 
from Eq. ()48|) that k_i^typ ~ i^\\,typY^^- I^ ^.Iso follows from our consideration above that 
typical values of k\\ are of order |g|. Combining the pieces, we find that the integral diverges 
as |g| ' , or 

6x{q)cxT\n)\q\'/^ (50) 

Similarly, at finite T we obtain 

6x{T)cxT\n)\T\'/' (51) 

And the specific heat 

5C(T)/T ocr2(7r)|T|2/3 (52) 

For special inflection points along symmetry direction, the analogous consideration shows 
the angle integral yields l/(fc±,typ)^- At the same time, it follows from Eq. (j^ that k±^typ ~ 
i^\\,typY^'^ ~ \q\^^'^- Then the angle integral diverges as \q\^^'^, and 

5x(g)ocr2(7r)|g|i/2 (53) 

while 

Sx{T)^T\n)\T\'/' (54) 

and 

5C(T)/T ocr2(7r)|T|i/2 (55) 

The results, Eqs. (j3Ul - (fS^ ). differ from the results by Fratini and Guinea^. They obtained 
x(T) oc TlogT for a generic inflection point, and xiT) oc T^^^logT for a special inflection 
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point. In our calculations, a similar Tlog^T behavior for a generic case would hold if the 
non-analytic terms were coming from vertices T{6) with arbitrary 6 rather than 6 = n. Then, 
e.g., the coefficient for the |g| term in the spin susceptibility would be given by a double 
integral J d6i\ko{6i)\ J (i^2|^o(^2)|- Each integral diverges logarithmically, and the correction 
would then scale as Tlog^T. The emergence of the anomalous power of temperature or 
momenta for a generic inflection point in our calculations is the direct consequence of the 
one-dimensionality of the relevant interaction. We see therefore that the anisotropy of the 
Fermi surface is an ideal tool to probe the fundamental ID nature of the non-analyticities 
in a Fermi liquid. 

VII. APPLICATION TO Sr2RuOi 

A crucial and so far unresolved question related to the results reported here and in pre- 
vious papers is the observability of the effects. Evidence for the T^lnT nonanalyticities 
expected in three dimensional materials have been observed in the specific heat of ^He^ 
(indeed this observation played a crucial role in stimulating the theoretical literature) and 
similar effects have been noted in the specific heat of the heavy fermion material UPt^. 
More recently, a linear in T behavior has been observed in the specific heat of fiuid mono- 
layers He^ adsorbed on graphite^'^, but to our knowledge no evidence for nonanalytic terms 
in the susceptibility has been reported. 

We consider here Sr2Ru04, a highly anisotropic layered compound for which detailed 
information about the shape of the Fermi surface, the quasiparticle mass enhancements, 
the susceptibility and optical conductivity is available2*i^^. These data imply^ that the 
material is "strongly correlated", in the sense that Fermi velocities and susceptibilities are 
substantially renormalized from the predictions of band theory. The data also suggest that 
the dynamical self energy is at most weakly momentum-dependent, because the shapes of the 
Fermi surface deviate only slightly from those found in band structure calculations, implying 
that the self energy has a much stronger frequency dependence than momentum dependence. 
We extract Fermi surface shape and Fermi velocities from quantum oscillation data^ and 
cross-check with photoemission data where available. We use available susceptibility^ and 
optical data-^-' to estimate the interaction functions. 

It is also useful to consider the leading low-T behavior of the specific heat coefficient. 
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which in a Fermi hquid is given in terms of fundamental constants and a sum over bands of 
the average of the inverse of the Fermi velocity as 



T 3 

A 



(57) 



with 



47r2 \vF,xi 

For Galilean-invariant fermions, it is customary to use the relation (we use units where 
h = 0) vp = kp/m to define a band mass 

mA = 2Tih (58) 

so C /T = (7r/3)^;^mA. To obtain the specific heat in conventional units {mJ/mol/K^) 
one must multiply by k'^ and by the Avogardo number. 

We begin with the results for the Fermi surface shape and Fermi velocities. In Sr2Ru04, 
the relevant electrons are the three t2g symmetry Ru d-orbitals and there are accordingly 
three bands at the Fermi surface, conventionally labeled as a, (3, 7. The Fermi surface shape, 
shown in Figgis well described by the two dimensional tight binding model 

e-yikx, ky) = — 2ti^ {cos{kx) + cos(ky)) — 4^27 {cos{kx)cos{ky)) — Eq^ (59) 

£a,l3 = -(tla + ha) {cOS^k^) + COs{ky)) - £0x2 

±Y ((^la - t2a){cos{kx) - cos{ky)f + lQtl^sin{kxY sin{kyY (60) 

(couplings in the third dimension are an order of magnitude smaller). 

A detailed quantum oscillation study has been performed by Bergemann and 
collaborators^. These authors present in Table 4 of their work a tight-binding parametriza- 
tion which reproduces the shape of the Fermi surface. They also present results for the 
mass enhancements in each Fermi surface sheet, which may be converted into experimental 
estimates for I\. The shape, of course, does not depend on the magnitudes of the tight bind- 
ing parameters. We accordingly rescale these in order to obtain velocities (more precisely, 
integrals I\) corresponding to the data reported by Bergemann et al^. 

A few remarks about the velocities and masses are in order. First, the calculated 7 band 
properties depend very sensitively on how close the 7 (a;?/)-derived band approaches the 
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TABLE I: Tight binding band parameters (in [eV]) which reproduce the shape and, approximately, 
the Fermi velocities of the three bands at the Fermi surface of Sr2RuOi. Parameters are taken from 
Table 4 of Refi and then renormalized to produce sheet-dependent quasiparticle mass enhancements 
approximately consistent with experiment. Last column: mass parameter computed using Eqs ()57() . 
(EOl). 



band 


£0 


ii 


t2 


H 




a 


0.13 


0.13 


0.013 0.02 2.5 


a 


0.16 


0.15 


0.013 0.02 5.8 


7 


0.012 0.079 0.032 
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van Hove points (7r,0), (0,7r). Published band calculationsiSiiLiS show wide variations in 
the position of the the singularity relative to the Fermi level. Second, the 'mass' derived 
from the specific heat involves both the velocity and the geometrical properties of the Fermi 
surface. The mass for the a band is small because of its small size, even though its velocity 
is relatively small. Third, and most important, the curvature of the a,/3 bands depends 
very sensitively on the parameters t^a.Ha] the velocities also depend somewhat on these 
parameters. A recent angle-resolved photoemission experimenlji^ reports that the a band 
Fermi velocity at the zone face crossing point is Va = 1.02eV^ — A; the parameterization used 
here gives an essentially identical value. 

We now turn to the Landau interaction function. A complete experimental determina- 
tion is not available, but considerable partial information exists. Bergemann and co-workers^ 
have determined, for each band, the spin polarization induced by a uniform external mag- 
netic field, so the "L = 0" spin channel Landau parameters may be estimated. Optical 
conductivity data^^ provide some information on the "L = 1" spin-symmetric channel cur- 
rent response. General arguments suggest that the charge compressibility is only weakly 
renormalized in correlated oxide materials, allowing a rough estimate of the "L = 0" charge 
channel interaction. We will use this information to estimate the scattering amplitudes and 
hence the nonanalytic terms in the susceptibilities. These estimations are certainly subject 
to large uncertainties, but we hope they will give a reasonable idea of the magnitude of the 
effects. 
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In the three-band material of present interest the interaction function is a symmetric 
3x3 matrix with components F'^''' labelled by orbital or band indices, which should then 
be decomposed into charge (symmetric) and spin (antisymmmetric) components and into 
the angular harmonics appropriate to the tetragonal symmmetry of the material. We begin 
with the isotropic "L = 0" spin channel. We assume (consistent with the usual practice in 
transition metal oxides) that the deviations from 0(3) symmetry, while crucial for electronic 
properties such as the band structure and conductivity, are not crucial for the local interac- 
tions, which arise from the physics of the spatially well localized d electrons. This implies 
that the interactions are invariant under permutations of orbitals, so that it is reasonable 
to assume that the two-particle irreducible spin channel interaction takes the simple Slater- 
Kanamori form with two parameters, which we write as F"'" = U^ff and F"^** = Jeff- Thus 
the physical static susceptibilities are given by 



-1 



X={Xo'+nUeff,Jeff)y (61) 

and fix the parameters U^f f and J^ff by comparing measured susceptibihties to the values 
predicted by the renormalized tight binding parameters. 

Ref^ presents (as mass enhancements) data for the spin susceptibility of each band (ob- 
tained from the spin splitting of the Fermi surfaces); finding x"'"/Xo'" ~ 1-2, X^'^ /Xo ~ 1-3 
and x?''^ IXq'^ ~ 1-6- We estimate Ueff ~ 0.033 and Je// = —0.008, where xo is the suscep- 
tibility implied by the quantum oscillations Fermi surface and mass. This implies that the 
dimensionless Landau interaction parameters (in the limit ujjk — > limit A°^ = T°^^^^/x^(^ 

^-0.053 0.040 0.10 ^ 

0.04 -0.13 0.17 (62) 

\ 0.10 0.17 -0.40 , 

The uncertainties in the off diagonal components are large, perhaps 50%, but because 
the interactions enter squared, the contribution of the off diagonal components is not very 
significant. The dominant term is the 7 — 7 band interaction, as expected because it has 
the largest mass and the largest susceptibility enhancement, but that all of the other con- 
tributions taken together make a non-negligible contribution to the interaction. Finally, we 
note that the spin channel renormalizations are not large, so use of the second order result 
is not unreasonable. 
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spin 



We now turn to the charge channel, beginning with the compressibihty. There is no 
experimental information available. However, it is generally believed that for systems, such 
as transition metal oxides, with strong local interactions the total charge susceptibility is 
not strongly renormalized, so that the Landau parameter acts to undo the effects of the 
mass enhancement. Further, if the J (orbital non-diagonal) component of the interaction is 
not too small relative to the U (orbital diagonal component) then a residual interaction acts 
to shift the levels such that the ratio of occupancies of each of the three t2g orbitals remains 
constant under chemical potential shifts. Taking as unrenormalized value the susceptibilities 
folUowing from the tight binding parameters given in Re^ we then obtain 

^ 0.40 0.034 0.026^ 
0.034 0.8 0.026 



A 



Sfl 



(63) 



\ 0.026 0.026 0.80 / 



with again considerable uncertainty in the off diagonal components. 

Finally, we turn to the current renormalization. The optical conductivity is commonly 
presented in the extended Drude form 



I.I?, 



band 



= ru^Zpi (64) 

where c is the mean interplane spacing and vn* jvn has the meaning of an optical mass 
enhancement defined with respect to a reference value determined by Dfjand- In a Fermi 
liquid at low temperatures, T{VL -^ 0) is very small and (assuming tetragonal symmetry) 



-^^^^, -^O-^j §MM ^ H- (m)'^ ,...,.)| (l . f ) ,65, 

A ' 

Note that the numerical value of the mass enhancement m* /mP depends on the choice of 
reference value Diyand but that D is a physically meaningful quantity determined directly 
from the data. 

The room temperature conductivity of Sr2Ru04 has been measuredJ^. These authors 
chose the value e'^Dhand/^c (which they denote as Up/An) to correspond to cijp ^ 8 x lO^cm"^ 
and find m* /m{Q -^ 0) (which they denote as A) to be ^ 3.5. This implies that D ^ 0.13eV^, 
somewhat smaller than the value 0.18 obtained from EqlHHl implying that the average over 
all bands is F^^ ^ —0.55. The temperature dependence of D in Sr2Ru04 has not been 
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measured, but it seems reasonable that D should decrease as T decreases, implying a further 
increase in the magnitude of F^^ . Determining the temperature dependence of the optical 
mass is therefore an important issue. 

Re^ presents, as masses, data for the cyclotron resonance frequencies for the different 
Fermi surface sheets. These masses should be essentially equivalent to the D values quoted 
above. Bergemann et. al. emphasize that the frequencies are subject to large errors, and 
that the results should be regarded as tentative. The quoted cyclotron masses correspond to 
D values about a factor of two larger than those implied by the measured fermi velocities, 
and about a factor of 3 larger than the values inferred from the optical data. In view of 
the stated large uncertainties in the measurement and the qualitative incosistency with the 
optical data, we disregard the cyclotron resonance measurements here. 

Now, the crucial object for the specific heat is the backscattering amplitude. A negative 
F^^ implies a positive backscattering amplitude, so as a rough approximation to the effects 
of the current channel Landau renormalization we add the interaction corresponding to 
F^^ = —0.6 to the diagonal components of EqlH!^ 

The crucial points emerging from this analysis are that the reducible interactions in the 
charge channel are of order unity, whereas those in the spin channel are somewhat smaller, 
implying a larger nonanalyticity in the specific heat than in the susceptibility. Substituting 
the interaction amplitudes into Eqs. |2Zl ^S and performing the fermi surface averages then 
yields the following estimates 

7(T) = 36mJ/mol-K^{l-0mi5T[K]) (66) 

x{T)[St /Volume] = 1.5 x 10"^ (1 - .00001T[ii']) (67) 

The small magnitude of the corrections (especially to the spin susceptibiltiy) follows from 
the small prefactors in Eq (PT|) and the not too large Landau renormalizations. The size of 
the effect is increased by the relatively small curvatures of the a and, especially, (3 bands, and 
we note that substantial increases in the coefficients occur if the mixing coefficients t2a, tsa 
in EqEOlare reduced. We expect the results to be valid above a (still not well determined) 
scale probably ~ 1 — 2K at which the Fermi surface warping becomes important enough 
to make the material three dimensional, and below the scale at which Fermi liquid theory 
breaks down, and we see that temperatures of order lOK lead to 20% deviations in the value 
of the specific heat coefficient and to 1% changes in x- 
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Replacing Sr by Ca leads to a dramatic (factor ~ 100 in SrQjjCaijjRu04) enhancement 
of the susceptibility. It seems likely that this increase is not due to a decrease in the fermi 
velocities, but must be interpeted as a dramatic increase in the spin Landau parameter, 
suggesting perhaps that nonanalytic T-dependence of x might be more easily observed in 
Ca doped materials, although in this case disorder effects would need to be considered. 

VIII. CONCLUSIONS 

In this paper, we studied non-analytic terms in the spin susceptibility and specific heat 
in 2D systems with anisotropic, non-circular Fermi surfaces. For systems with circular 
Fermi surfaces, the non-analytic terms in Xs(g, T) and C{T)/T are linear in max{q,T). 
We argued that the anisotropy of the Fermi surface serves as a testing ground to verify 
the theoretical prediction that the non-analytic terms originate from a single ID scattering 
amplitude which combines two ID interaction processes for particles at the Fermi surface in 
which the transferred momenta are either or 2kF, and, simultaneously, the total moment 
is zero. We obtained explicit expressions for the non-analytic momentum and temperature 
dependences of the spin susceptibility and the specific heat in systems with non-circular 
Fermi surfaces and demonstrated that for the Fermi surfaces with infiection points, the the 
non-analytic temperature and momentum dependences are Xs oc max{q^/^^ T"^^"^)-, C{T)/T oc 
T^/^ in a generic case, and as Xs oc max(g^/^,T^/^), C{T)/T oc T^^"^ for the special cases 
when the infiection points are located along symmetry axis for the quasiparticle dispersion. 
We estimated the order of magnitude of the effects in the quasi two dimensional material 
Sr2Ru04^. 
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useful conversations. The research is supported by NSF Grant No. DMR 0240238 (AVC) 
and DMR 0431350 (AJM), and AJM thanks the DPMC at the University of Geneva for 
hospitality while this work was completed. 
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APPENDIX A: THE DETAILS OF THE EVALUATION OF EO. (|5^ 



For simplicity, we neglect band index, i.e., set v^ = v\ = Vk, and fco(^) — ^o(^) — ^oiM). 
Using (|M|). Eq. (^^ is re-expressed as 



8xLw{q) 



t\'k) fddkHk) 



647r" 



2txvI 



A /.A 



dxdy 



A J -A 



[{xy — e)^ log {xy — e)^ + {xy + e)^ log {xy + e)^ — 2x^y^ log x^y^] 



(Al) 

where e = (v^ ■ q)kQ{k)/vk- Rescaling x = \/\e\x, y = \/]e\y, substituting into (jAlj) and 
dropping irrelevant terms confined to high energies, we obtain 



SxLwiq) 



T\n) rd9kko{k)\e\ 



64v]-5 



2nvl 



(A2) 



where 



/A pl\ 
/ dxdy \{xy — tf log {xy — ef + {xy + ef log {xy + tf — 2{xyY log (xy)^] 



(A3) 



Introducing further x = v2rcos0/2 and y = v2rsin0/2, we rewrite Z as 

/»7r /"A 

Z = 2 d'ip dr ^{r sin — 1)^ log (r sin — 1)^ + (r sin + 1)^ log (r sin — 1)^ 

Jo Jo 

-2r^ sin^ log r^ sin^ 0] ( A4) 



Subtracting the irrelevant large r contribution 6 + logr^sin^0 from the integrand in ()A4|1 . 
we the universal part of Z in the form 

/ 1 \^ f'^ + V^^ ^ 

r^ sin^ log I 1 ^ — 1 + 2r sin log ' 



Z = 2 di) dr 

'0 Jo 



r^ sin^ - 



1- 



]_ 

r sin 4> , 



—6 + log 



(r^ sin^ - 1)2 
r^ sin^ 



(A5) 



One can make sure that the integral over r vanishes if we set the upper limit at A = oo. as 
the integrant depends on r only via r sin 0, the finite contribution to the integral comes from 
Asin0 = 0(1), i.e. from a narrow range of either near zero or near tt. The contributions 
from these two regions are equal. Resticting with the contribution from small 0, expanding 
sin 0^0 and introducing z = rcj) and t = A0, we obtain from ()A5|) 



Z = A 



dt 



dz 



z^ log 1 



1 



+ 2z log 



i + i 
1-i 



— 6 + log 



z^-l 



(A6) 
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Changing the order of the integration, we obtain for the universal part of Z 

2;2log(l-^j +22;log(^— fj - 6 + log 



-4 / dzlogz 
'o 



Z2 



Att 



2 



(AT) 



3 

Substituting this into ()A2|) and using the definition of e, we reproduce (jH^ . 

APPENDIX B: AN ALTERNATIVE EVALUATION OF 5xL^(g) 

In this Appendix we present a complementary evaluation of Sxiwil) using a somewhat 
different computational procedure. We again restrict to one band. The point of departure 
are Eqs. ()28|) and ()29|). which we re- write at T = as 



5xw(g) = -4/ / / / ^'^f^f^" r^n(/,(^)Go(k,.;)x 




(27r)6 
Go(k + I, cu + fi) Go(k + q + 1, cu + fi) Go(k + q, u) (Bl) 

We first integrate over internal momenta k and frequency a; in the fermionic propagator. 
Expanding the result in g^, we obtain 

SxLwiq) oc rV Hdnn f dH [ do,^^ — — ^n(/,n) (B2) 

Jo J J vf{9i) {lil - V fI cos if 

where 9i is the angle between 1 and k. Directing l^ and ly along and transverse to k and 
substituting the polarization operator we obtain 

1 



5xw(g)ocrV r d^^l^ I dU I dOi^ 
Jo J J vf{9i 



) {zn-VF{ei%r 



X 



di /■rf^M^i±^ i rB3i 

V V f{9i + 9) iQ - V f{9i + 9) {l^ COS 9 + ly sin 9) ^ ' 

[9 is the angle between two internal momenta p and k]. We now integrate over ly and then 
over 9. The full result for this 2D integral depends on particular forms of ko{9) and vf{9). 
However, we only need from the integral over dlyd9 the term which is non-analytic in the 
lower half-plane of Ix (this will allow us to avoid a degenerate pole at vfIx = i^) ■ One can 
easily verify that the non-analyticity comes from the integration near 9 = it which yields, 
instead of the second line in ()B3|) 

' 2vll9i~-l) ^°^ ['^ + ""^^^^ " ^^^^'l ^^^^ 
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For the Fermi surfaces with inversion symmetry (which we will only consider) ko{6i — tt) = 
ko{9i) and vf{Oi — tt) = vf{Oi) (we recall that vf{0) is the modulus of the Fermi velocity 
at a particular 9). Substituting this result into ()B3|) and extending the integral over l^ onto 
the lower half-plane, we obtain Eqn (j35p . 

We also verified that the same result could be obtained by evaluating the singular part 
of n(/, Q) by explicitly expanding near p = — k and expanding the dispersion Sp = S-k+i 
to second order in /. In this computation, one power of A;o(^) comes from expanding the 
dispersion, while the other comes from the Jacobean of the transformation from (Pk to 
dSkdO. 

APPENDIX C: REEVALUATION OF 5xLw[q) FOR AN ISOTROPIC FERMI 
SURFACE 

In this Appendix, we reconsider a previously published^ evaluation of 5xlw{i)- Although 
this evaluation leads to results identical to those we presented in the body of the paper 
for a circular Fermi surface, it apparently does not invoke the curvature explicitly. Here 
we deconstruct this analysis, showing how the curvature actually enters even when Fermi 
surface is circular. 

We begin from Eq ()B1|) . The analysis presented in the main body of the paper involves 
choosing a direction for I, and then performing the integral over k, which picked out points 
with a definite relationship to 1 and involved the curvature in a direct way, and finally 
integrating over /. On the other hand, the "conventional" analysis involves first fixing the 
direction of k, integrating over the magnitude of k and over g, and then averaging over the 
direction of k. In this method one expands e^, eu+i, e^+q and Cfc+z+g in ()B1|1 to linear order in 
the deviations from the Fermi surface as e^ = vpik — kp), e^+i = e^ + vplx, etc. Because the 
Green functions have been linearized the curvature apparently does not enter, in contrast 
to the previous derivation, where the dependence of the Green function lines on curvature 
was essential. 

Integrating over k and over the corresponding Matsubara frequency, and expanding the 
result in powers of q, we obtain at T = 0, neglecting regular terms 



6xLw{q) oc r\^ J dQQ J d/,— — 1-— J dlyU{l, Q) 



(CI) 
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The key point is that the curvature dependence is hidden in the polarizibihty 11, but in the 
circular Fermi surface hmit this dependence is hidden. To make the curvature dependence 
manifest we use the fact that only backscattering contributes and evaluate the polarization 
bubble n(/, Q) = J dPtduj'GQ{t, uj')Gq{1 + t, fi + uj') by expanding near t = — k. Introducing 
t + k = p and assuming that p is small, we expand the dispersions e^ = e_fc^+p and tt+i = 
e-fc^+p+i to second order in p: 

f P^ \ ( (P +1 ?\ 

et = e_fc,+p = -vf [P^ + ^j'^ ^w = e-fe,+p+i = -vf {px + lx+ ^ ^^ ^ 1 (C2) 

Substituting this expansion into the bubble and integrating over py we obtain 

Till n\ ■ ^^0 1 Aly-{in + VFlx) .^„^ 

n(/,n) = 1 log — 7^ — — (C3) 

^ 27r^Vply -Aly - {iVt + vfIx) 

where k^A/vF ~ kF is the upper limit of the integral over py. Integrating next n(/, Vt) over 
ly, we find the same branch cut singularity as in "conventional" approach 

// o 
diyiia, n) = ^ iog[zf] + vfU (C4) 

Tcvi 



Substituting this result into (jBlll and using the fact that (fk in (jBlll can be re- expressed 
as {kQ/vF]dekd9, we reproduce Eq. (jH^j) for a circular Fermi surface, and also reproduce Eq. 
(4.18) in ^5a], but with ko/vF instead of m. 

For completeness, we also show that Sx2kp{(l) in systems with a circular Fermi surface 
can also be obtained with and without the curvature. A "conventional" computation [5|a] 
expresses Sx2kp (q) in terms of the curvature. An alternative computational scheme involves 
the same "triad" method that we used in the main text. In this scheme, the original 
expansion near 2kF momentum transfer contains the curvature, but it disappears from the 
answer at the latest stage. Performing the same integrations over e^, the corresponding 
frequency and ly as in the main text, we find (keeping F = r(6')) 



//•oo 
deT\9) / dnn 
J Vp\n\ 



-VfW 

— 1 -0 T (C5) 

{i^ - iny (in -vf{Ix cose + ^8111^9))^^ ^ ' 

For cos 9 < 0, the two double poles are in different half-planes of Ix- Integrating over l^, we 
then obtain 

Sx2kA<l)^Vh d9T\9) / ^^,„ . — - (C6) 

Jn/2 Jvp\q\ {^2^ Sm^ 9 + lil{l - COS 9)^''^ 
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Since relevant Q ~ '^^fIo'I, the 6 integral is confined to ^ = vr. Expanding near vr we obtain 



l/2(^sin2^ + ifi(l-COSe))5/2 ' Vo (^x2-2zfi)5/2 Sfi^^^g 

Substituting this into ()C6|) . we find that k^ is canceled out, and 






Restoring the prefactor, we reproduce the same result as in the main text, but with mvp 
instead of /cq- 
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